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Abstract 

In this note, we discuss the basic elements that should appear in a gravita- 
tional system dual to a confining gauge theory displaying color superconductivity 
at large baryon density. We consider a simple system with these minimal ele- 
ments, and show that for a range of parameters, the phase structure of this model 
as a function of temperature and baryon chemical potential exhibits phases that 
can be identified with confined, deconfined, and color superconducting phases 
in the dual field theory. We find that the critical temperature at which the 
superconducting phase disappears is remarkably small (relative to the chemical 
potential). This small number arises from the dynamics, and is unrelated to any 
small parameter in the model that we study. We discuss similar models which 
exhibit flavor superconductivity. 
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1 Introduction 



Background 

Quantum Chromodynamics is believed to display a rich phase structure at finite tem- 
perature and chemical potential, with phase transitions associated with deconfinement, 
nuclear matter condensation, the breaking of (approximate) fiavor symmetries (which 
are exact in generalizations with equal quark masses and/or massless quarks), and the 
onset at high density of quark matter phases displaying color superconductivity (for 
reviews see for example [HEllSlllj). However, apart from the regimes of asymptotically 
large temperature or chemical potential, a direct analytic study of the thermodynamic 
properties of the theory is not possible. 

Even using numerical simulations, only the physics at zero chemical potential is 
currently accessible, since at finite /i the Euclidean action becomes complex, and the 
resulting oscillatory path integral cannot reliably be simulated using standard Monte- 
Carlo techniques. Current proposals for the phase diagram of QCD and related theories 
are largely based on qualitative arguments and phenomenological models. While these 
provide a plausible picture, it is possible that they miss important features of the 
physics. It would certainly be satisfying to have examples of theories similar to QCD 
in which the full phase diagram could be explored directly. 

The holographic approach 

A modern route to understanding properties of strongly coupled gauge theories, that 
would be otherwise inaccessible, is via the AdS/CFT correspondence, or gauge the- 
ory / gravity duality. This suggests that certain quantum field theories (usually called 
"holographic theories"), generally with large-rank gauge groups, are equivalent to grav- 
itational systems. By this correspondence, calculations of physical observables in the 
field theory are mapped to gravitational calculations; in many cases difficult strongly- 
coupled quantum mechanical calculations in the field theory (such as those required to 
understand the thermodynamic properties of QCD) are mapped to relatively simple 
classical gravity calculations. Optimistically, it may then be possible to find a theory 
qualitatively similar to QCD for which the physics at arbitrary temperature and chem- 
ical potential can be understood exactly via simple calculations in a dual gravitational 
system. 

By now, there are well-known examples in gauge-theory / gravity duality for which 
the field theory shares many of the qualitative features of QCD (see, for example [S]). 
Further, many of these theories have been studied at finite temperature and chemical 
potential, revealing phase transitions associated with deconfinement, chiral symmetry 
breaking, meson melting, and the condensation of nuclear matter. However, to date, 
most of the theories that can be studied reliably using dual gravity calculations have 
the restriction that the number of fiavors is kept fixed in the large limit. In such 
theories, the physics at large chemical potentials is known to be qualitatively different 
than in real QCD. For example, at asymptotically large chemical potential, theories 
with large Nc and fixed Nf are believed to exhibit an inhomogeneous "chiral density 
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wave" behavior [6l [7], rather than the homogenous quark matter phases predicted 
for finite Nc and Nf. In order to find examples of holographic theories which most 
closely resemble real QCD at finite chemical potential, one should therefore attempt to 
find examples of calculable gravitational systems corresponding to theories with finite 
Nf/Nc- This situation presents some technical challenges, as we now review. 

In the well-known examples of holographic gauge theories, the addition of fiavor 
fields in the field theory corresponds to adding D-branes on the gravity side [8j . Quarks 
correspond to strings which have one endpoint on these D-branes, while mesons cor- 
respond to the quantized modes of open strings which begin and end on the branes. 
The configurations of these D-branes in theories with finite Nf and large Nc are deter- 
mined by finding action-minimizing configurations of the branes on a fixed background 
geometry. On the other hand, in order to have Nf of order Nc in a large Nc theory, we 
need a large number of these fiavor branes, and these will back-react on the spacetime 
geometry itself. For Nf ~ Nc, there are as many degrees of freedom in the fiavor fields 
as there are in the color fields (gauge fields and adjoints), so it is natural to expect that 
the back-reaction will be so significant that in the final description the fiavor branes 
themselves will be completely replaced by a modified geometry with fiuxes (in the same 
way that the branes whose low-energy excitations give rise to the adjoint degrees of 
freedom do not appear explicitly in the gravity dual description of the field theory). 

There has been significant progress in understanding the back-reaction of fiavor 
branes, with some fully-back reacted analytic solutions available (for a review see [9J), 
but so far, there has not been enough progress to fully explore the phase structure of a 
QCD-like theory with finite Nf/Nc- In particular, as far as we are aware, color super- 
conductivity phases have not been identified previously in holographic field theories]^ 

Quark matter from the bottom up 

In this paper, we aim to come up with a holographic system describing a confining 
gauge theory that does exhibit a quark-matter phase with color superconductivity at 
large chemical potential. However, motivated by recent condensed matter applications 
of gauge/gravity duality (see, for example [13j), we will avoid many of the technical 
challenges described above by taking what is known as a "bottom up" approach. Rather 
than working in a specific string theoretical model which takes into account the back- 
reaction of fiavor branes, we will make an ansatz for the ingredients necessary for such 
a model to describe the relevant physics. We study the simplest possible gravitational 
theory with this minimal set of features, with the hope that it captures the qualitative 
physics of interest. We will indeed find that even this simple theory exhibits many of 
the expected features. 

^However, see [12] for a possible manifestation of the related color-flavor locking phase in a holo- 
graphic system. 
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Ingredients 

We wish to construct a gravitational theory to provide a holographic description of a 
four-dimensional confining gauge theory on Minkowski space with Nf ~ flavors. On 
the gravity side, the Minkowski space will appear as the fixed boundary geometry of our 
spacetime, but we must have at least one extra dimension corresponding to the energy 
scale in the field theory. Since the field theory has a scale (the QCD or confinement 
scale), the asymptotic behavior of the solution must exhibit an additional scale relative 
to the asymptotically AdS geometries that appear in gravity duals of conformal field 
theories. In the simplest examples of gravity duals for confining gauge theories, this 
scale is provided by the size of an additional circular direction in the geometry]^ Thus, 
we will work with a gravitational system in six dimensions whose boundary geometry is 
R^'^ X S^. We will assume that the asymptotic geometry is locally Anti-de-Sitter space, 
so the confining gauge theory we consider arises from a five-dimensional conformal field 
theory compactified on a circle. When we study the theory at finite temperature, there 
will be an additional circle in the asymptotic (Euclidean) geometry, the Euclidean time 
direction whose period is 1/T. 

The gauge theories we are interested in have at least one other conserved current, 
corresponding to baryon (or quark) number. By the usual AdS/CFT dictionary, this 
operator corresponds on the gravity side to a U{1) gauge field in the bulk. The asymp- 
totic value of the time component for this gauge field corresponds to the chemical 
potential in our theory, while the asymptotic value of the radial electric flux corre- 
sponds to the baryon charge density in the field theory. For a given chemical potential, 
the minimum action solution will have some specific value for the flux, allowing us to 
relate density and chemical potential. 

The color superconductivity phases believed to exist at large density in QCD and 
related theories are usually characterized by condensates of the form (ipip), bilinear 
in the quark fields ip, which spontaneously break the U{N) gauge symmetry, and the 
U{1)b global symmetry. Naively, we would want to model such operators by a bulk 
charged scalar field corresponding to the condensate. However, bulk fields always 
correspond to gauge-invariant operators, while by definition the ipip bilinears which 
break the gauge symmetry are not gauge-invariant (in fact, there is no way to make 
a singlet from two fundamental fields, except in the case of SU{2)). Additionally, the 
simplest gauge-invariant operators charged under U{1)b involve ip fields and have 
dimension of order A^, thus our holographic dual theory should have no light scalar 
fields charged under the U{1)b gauge field. 

The correct way to understand the condensation of the ipip bilinears is as an example 
of spontaneously broken gauge symmetry (as in the Higgs mechanism), rather than as a 
phase transition characterized by some gauge-invariant order parameter. Nevertheless, 
the transition to color superconductivity can be characterized by the discontinuous 
behavior of gauge-invariant operators, which are of the form . Such operators 

are gauge invariant and neutral under the U{1)b, and therefore should correspond to 

■^There are other possibilities here, as we mention briefly in the discussion section. 
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an uncharged scalar field in the bulk with dimension of order l|^ 

Combining everything so far, we want to study gravity in six dimensions with 
negative cosmological constant and boundary geometry _R^'^ x with a U{1) gauge 
field and a neutral scalar field. The simplest action for this system i^ 

/ d'x^ |7^ + g - - - , (1) 

where we include one tunable parameter, the mass m of the scalar field, which de- 
termines the dimension of the corresponding operator in the dual field theory. More 
generally, we could consider other potentials for the scalar field, or a more complicated 
action (e.g. with a Chern-Simons term or of Born-Infeld type) for the gauge field, but 
we restrict here to this simplest possible modelF] 



Results 



Starting with the model ([T]), we have explored the phase structure by minimizing the 
gravitational action for specific values of temperature (corresponding to the asymp- 
totic size of the Euclidean 5*^ direction) and chemical potential (corresponding to the 
asymptotic value of Aq). Our results for the phase diagrams are shown in figures Ipp^ 



For small /i, we find a confined phase at low-temperature and a deconfined phase at 
high temperature, with the scalar field uncondensed in each case. However, increasing 
fj, at zero temperature, we find (setting Lacis = 1) for — ^ < < — 5 a transition 
to a phase with nonzero scalar condensate (on a geometry with horizon) and finite 
homogeneous quark density, as expected for a color superconductivity phase. Increas- 
ing the temperature from zero, we find a transition back to the deconfined phase at 
a remarkably low temperature; for example, at = —6, the critical temperature at 
which superconductivity disappears is 

T//i ~ .00006333 . 

The tendency for the scalar field to condense at low temperatures for the range of 
masses above can be understood in a simple way, as explained for example in [HI [15] . 
In d + 1-dimensional anti-de Sitter space with anti-de Sitter radius L, the minimum 

^As emphasized to us by Andreas Karch, a gauge invariant operator of the form O4 = tpip{ipipy 
can be written as a sum of terms OaOa where each Oa ~ (V'^^)q is gauge invariant (and a represents 
flavor/Lorentz indices). Thus, ©4 is something hke a double-trace operator. In a large TV theory, 
factorization of correlators implies that the expectation value of O4 can be calculated classically from 
the Oa expectation values (up to corrections). Thus, discontinuous behavior of O4 should be 
directly related to discontinuous behavior in the simpler gauge-invariant operators Oa (which also 
have no baryon charge), so it may be more appropriate to think of the scalar field in our model as 
being dual to one of these simpler operators. 

^Since we will also consider the case of a charged scalar field, we have written the action using 
standard normalizations for a complex scalar, but we will take the scalar to be real in the uncharged 
case. 

^For another approach to modeling the QCD phase diagram by an effective holographic approach, 
see for example [lOj [11] . 
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Figure 1: Phase diagram of our model gauge theory with = —6, R = 2/5. Region 
in dashed box is expanded in next figure. 

mass for a scalar field to avoid instability is m'^BF — — <i^/(4L^). The minimum action 
solution for large chemical potential in the absence of any scalar field is a planar 
Reissner-Nordstrom black hole solution with one of the isometry directions periodically 
identified. In the limit of zero temperature, the near horizon region of this black hole 
has geometry AdS'^ x R'^, with the radius of the AdS'^ equal to L2 = L/\/20. Thus, 
in the near-horizon region, there will be an instability toward condensation of the 
scalar field if < —l/lAL"^) = — We thus have a range (setting L = 1) of 
—25/4 < < —5 for which the scalar field tends to condense in the near-horizon 
region but is stable in the asymptotic region. Numerical simulations verify that we 
indeed have scalar field condensation for precisely this range of masses. 

While there is no guarantee that the gravitational system we study has a legitimate 
field theory dual, "top-down" gravitational systems corresponding to fully consistent 
field theories must have the same basic elements (usually with additional fields and a 
more complicated Lagrangian). The fact that the expected physics emerges even in 
our stripped-down version suggests that quark-matter phases will be found also in the 
complete models, once back-reaction effects are under control. Optimistically, qualita- 
tive features that we find in the bottom-up model (such as the extremely low transition 
temperature between superconducting and deconfined phases) may be present also in 
more complete holographic theories. In this case, our simple model may provide novel 
qualitative insights into fully consistent QCD-like theories. 
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Figure 2: Phase diagram of our model gauge theory with = —6, i? = 2/5. Region 
in dashed box is expanded in next figure. 




Figure 3: Phase diagram of our model gauge theory with = —6, R = 2/5. The 
dashed curve represents the phase boundary in theory without a scalar field. 
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Charged scalar 

While less relevant to color superconductivity, it is also interesting to explore the 
physics of our model when we make the scalar field charged under the gauge field. In 
this case, the scalar field corresponds to a gauge-invariant operator in the field theory 
that is charged under the U{1) associated with A, and the kinetic term for the scalar 
field is modified in the usual way as d^ip — )■ d^ip — iqAfj^ip. As we have argued above, this 
symmetry cannot be U{1)b, but could be another fiavor symmetry, such as isospin in 
a model with two or more fiavors. The fiavor superconductivity associated with meson 
condensation was studied previously in the holographic context (with finite Nf), for 
example in [HI |l7l [18]. Our results are qualitatively similar to the ones obtained in 
those studies, and we leave more detailed comparison for future work. 

In section 5 below, we determine the phase diagram for various values of q and 
m. The same system was studied for the 2+1 dimensional case in jTH] and origi- 
nally in [2U] for the case of large q. The application there was to holographic insula- 
tor/superconductor systems, but the intriguing resemblance of the phase diagrams in 
those papers to QCD phase diagrams partially motivated the present study. 

2 Basic setup 

In this paper, we consider holographic field theories with a conserved current J^, as- 
sumed to be a baryon current (or isospin current when we consider charged scalar 
fields) and some gauge-invariant operator O whose condensation indicates the onset of 
(color or fiavor) superconductivity. We would like to explore the phase structure of the 
theory for finite temperature T and chemical potential /x; that is, we would like to find 
the phase that minimizes the Gibbs free energy density g = e — Ts — fip, where e, s, 
and p are the energy density, entropy density, and charge density in the field theory. 
We can also ask about the values of e, s, p, and (O) as a function of temperature and 
chemical potential. 

As discussed in the introduction, our holographic theories are defined by a dual 
gravitational background which involves a metric, U{1) gauge field, and scalar field, 
with a simple action 

Jd'xV^g[n+f^-\F'-\d,^\'-m'\i.\'] . (2) 

We choose coordinates [t, x, y, z) for the non-compact field theory directions, w for 
the compact field theory direction, and r for the radial direction. We take boundary 
conditions for which the asymptotic (large r) behavior of the metric is 

ds^ (-dt^ + dx'^ + dy"^ + dz^ + c?w^) + dr^ , 

where w is taken to be periodic with period R. To study the theory at finite tempera- 
ture, we take the period of t = it in the Euclidean solution to be 1/T. 
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The equations of motion constrain the gauge field to behave asymptotically as 



A — _L 
■r^i, — Ctjy — + . . . . 

Since Ai, is assumed to be the field corresponding to the conserved baryon current 
operator J", in the field theory, the usual AdS/CFT dictionary tells us that is 
interpreted as the coefficient of the J'^ in the Lagrangian (i.e. an external source for 
the baryon current) while is interpreted as the expectation value of baryon current 
for the state corresponding to the particular solution we are looking at. To study 
the theory at finite chemical potential without any external source for the spatial 
components of the baryon current, we want to take 

ay = {jj, 0, 0, 0) . 

The scalar field equations of motion imply that asymptotically 

ip = — + — + ■■■ (3) 

where 

The holographic field theories we consider are defined by assuming ipi = 0. In this 
case, A+ gives the dimension of the operator dual to ip^ In this case, 1^2 (which will 
be different for solutions corresponding to different states of the field theory) gives us 
the expectation value of the operator O in the field theory. 

By the AdS/CFT correspondence, the field theory free energy corresponds to the 
Euclidean action of the solution. Thus, to investigate the field theory state which 
minimizes free-energy for given T and /i, we need to find the gravitational solution 
with boundary conditions given above which minimizes the Euclidean action. Note 
that we only consider solutions with translation invariance in t, x, y, z, and w. It 
would be interesting to investigate the possibility of inhomogeneous phases (or at least 
the stability of our solutions to inhomogeneous perturbations) but we leave this as a 
question for future work. 



Calculating the action 



In order to obtain finite results when calculating the gravitational action for a solu- 
tion, it is important to include boundary contributions to the action. In terms of the 
Lorentzian metric, gauge field and scalar, the fully regulated expression that we require 
is [131 



S 



lim 



r<rM 




DM' 



m 



^For a certain range scalar field masses in the range —(P/A < m? < —cP/A + 1, it is also consistent 
to define a theory by fixing ■02 = 0. In this case, the dimension of the dual operator is A_. We consider 
this case briefly in section 4.2. 
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+ / d'^xy/^ I -2K + 

Jr=rM [ 



where 



A_ = ^ - + 4m2 . 



Here, 7 is the metric induced on the boundary surface r = tm, and is defined as 

K = r^^^^n, , 

where is the outward unit normal vector at r = Tm- The scalar counterterm here is 
the appropriate one assuming that our boundary condition is to fix the coefficient of 
the leading term in the large r expansion of ip. Since we are setting this term to zero, 
it turns out that the counterterm vanishes in the tm — >■ 00 limit. 
For all cases we consider, the metric takes the form 



ds^ = j^^^i + 9oo{r)dt'^ + grr{r)dr^ + gww{r)dv? . 



(4) 



Assuming the Einstein equations are satisfied, we can show (by subtracting a term 
proportional to the xx component of the equation of motion) that the integrand in the 
first term may be written as a total derivative with respect to r 



-g{n + 



Using 
we have 



d{d-l) 
L2 



m 



n,. 



(0, ...,0, 



dr 



J9r 



K = r^^^n. 



= 7'''^{-r>.} 

7"" 



1 ugav 
9 



2" dr 
1 „..57mv 



-7 



1 (91n(^/=^) 



/9r 



dr 



so that 



2 9v/=7 
faZ dr 



Our final expression for the action density is 

p I'M 

S/Vd^ —^9 



rgr 



2 ay=7 2{d-l) , ^ 

+ T V -7 



'9r 



dr 



r=rM 



(5) 
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Action in terms of asymptotic fields 

It is convenient to rewrite the expression ([S]), in terms of the asymptotic expansion of 
the fields. For the ansatz (|4]), and the boundary conditions appropriate to our case, we 
find 

(3) 

2 , 9tt I 

9tt = -r + — + ... , 

_ 1 ^ 

9rr n ~^ 7 + • • • , 

o(3) 

„ = + + 
yww I ' 3 ' 

^(3) 

^ = 7r + --- ' 

^ = ^-3r^ + --- • 

Inserting these expansions into our expression above for the action we find that (as- 
suming the term at r = ro vanishes) 

S = 5g^^l + Agil'>-5gl? . 

However, using the equations of motion, we find that g'l'^l + g^^) — g^^ = 0, so we can 
simplify to: 

S = -gil^ . (6) 

Numerically, it can be a bit tricky to read off gll^ because there is also a term in 
the expansion of grr- But using the equations of motion, we can find 

From this, it follows that the combination 

3 

-r'^grrir) + - ^{7 + m^)r^il?{r) 

behaves like 

+ 0{l/T^) . 

So, we can numerically evaluate the action by taking 

S ^ -rlgrr{r^) + - ^(7 + m2)r^?/'^(rj , 
where r* is taken to be large but not too close to the cutoff value. 

3 Review: i/^ = solutions 

We begin by considering the solutions for which the scalar field is set to zero. 
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3.1 AdS Soliton solution 



At zero temperature and chemical potential, the simplest solution with our boundary 
conditions is pure AdS with periodically identified w. However, assuming antiperiodic 
boundary conditions for any fermions around the w circle, there is another solution 
with lower action. This is the AdS soliton [23], described by the metric (setting L = 1) 

ds^ = f-rft^ + dx'^ + dy'^ + dz^ + fir) dw'^) + , , (7) 



where 



7^5 



As long as we choose the period 2nR for w such that 

the solution smoothly caps off at r = tq. This IR end of the spacetime corresponds 
in the field theory to the fact that we have a confined phase with a mass gap. The 
fiuctuation spectrum about this solution corresponds to a discrete spectrum of glueball 
states in the field theory. 

Starting from this solution, we can obtain a solution valid for any temperature and 
chemical potential, by periodically identifying the Euclidean time direction and setting 
Aq = n everywhere. Using (|6| we find that the action for this solution is 



Ssol — '"o — 



5R 



The negative value indicates that this solution is preferred over the pure AdS solution 
with action zero. 



3.2 Reissner-Nordstrom Black Hole Solution 

For sufficiently large temperature and/or chemical potential, the AdS soliton is no 
longer the 1^ = solution with minimum action. The preferred solution is the planar 
Reissner-Nordstrom black hole, with metric 

ds'^ = (-dff{r) + dx^ + dy^ + dz^ + dw"^) + , (10) 

V ^ T f i T] 



where 



the scalar potential is 



0(r) = fiil- 
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Figure 4: Phase diagram without scalar field, in units where R = 2/5. 

and w is periodically identified as before. 

This solution has a horizon at r = r+. The temperature of the solution (determined 
as the inverse period of the Euclidean time for which the Euclidean solution is smooth) 
is given in terms of r_|_ by 

T = 1. Ur, - . ,12) 

47r V + 8r+; ^ ^ 

From ([6]), we find that the action for this solution is 



Srn = -rl 1 + 



Thus, we find that the black hole solution has lower action than the soliton for 




where r+ is determined in terms of T and /i by (12). This defines a curve in the T — ^ 
plane that begins on the = axis at T = l/(27ri?) and curves down to the T = 
axis at = 2^^/^° /{5^/'^3^/^R) ^ 4.3547/(27ri?), as shown in figure g 

As usual, the existence of a horizon in this solution indicates that the corresponding 
field theory state is in a deconfined phase |26j . 

In the next sections, we consider solutions with nonzero scalar field. We will find 
that for large fi there exist solutions with nonzero scalar field that have lower action 
than the solutions we have considered, so the phase diagram of figure |4] will be modified. 
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4 Neutral scalar field: color superconductivity 



In the case of a neutral scalar field, our simple model has no explicit source for the 
gauge field in the bulk, so homogeneous solutions with a non-trivial static electric field 
(corresponding to a non-zero baryon number density in the field theory) necessarily 
have a horizon from which the flux can emerg^ 

To look for solutions of this form, we consider the ansat^ 



At 



dr 



~g{r)e-''^'^dt^ + + r\dw^ + dx^ + dy^ + dz' 



The scalar and Maxwell's equations that follow from the action ^ are 

'4 x' 




9 
■ 



-V^ = 



while the Einstein equations are satisfied if 





5r = 



These have two symmetries: 



^(r) = i^iar) , 0(r) = -0(ar) , xir) = xiar) 

a 

arising from the underlying conformal invariance, and 

X = X + A , . 



9[r) 



-^g[ar 



_ A 

e 2 ( 



(13) 
(14) 

(15) 
(16) 

(17) 



We would like to find solutions with a horizon at some r = r+. The electric potential 
must also vanish at the horizon, and we are looking for solutions for which the leading 



falloff t/'i in rt3| vanishes for the scalar. Also, multiplying the first equation (13) g 



^In a more complete model, the source might be provided by some non-perturbative degrees of 
freedom in the theory, such as the wrapped D-branes that give rise to baryons in the Sakai-Sugimoto 
model. 

^We could have considered a more complicated ansatz, with an extra undetermined function in 
front of d'up' . However, it is plausible that as for the = solution, the minimum action solution for 
the case where the w circle does not contract in the bulk is a periodic identification of the solution 
with non compact w and rotational invariance in the x, y, z, w directions. 
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and evaluating at r = r+ fixes in terms of ipir.^) and g'{rj^). Altogetlier, our 

boundary conditions are 

^?(r+) = 0, 0(r+) = O, X(oo)=0, = , 

and 



40r+ - 2m?r\ip'^{r+) - r+e>^(''+)(0'(r+))2 ' 

The remaining freedom to choose r+ and leads to a family of solutions with 

different T and /i. Explicitly, we have 

Note that solutions with the same T/^ are simply related by the scaling symmetry 



(17) 



4.1 Numerical evaluation of solutions 



To find solutions in practice, we can make use of the symmetries (17) to initially set 



r_|_ = 1 and x(0) = and solve the equations with boundary conditions 

(7(1) =0, x(0) = 0, 0(1) = 0, <p\l) = E,, V'(l)=^o 

and 

^'(1) 



40 - 2m2^g - El 



We can integrate the and x equations explicitly to obtain 

x(r) = — f df-f 



leaving the remaining equations 



2 \dr ^ 
Ji 



r A g g 

r 4 8r' 4 

We use Eq as a shooting parameter to enforce t/'i = 0, and find one solution for each 



ipQ. From these solutions, we apply the symmetry (18) with A = — x(oo) to restore 



x(cxd) = and finally use the symmetry (17) to scale to the desired temperature or 
chemical potential. 

Using this method, we find that solutions exist for scalar mass in the range —25/4 < 
< 5, which is exactly the range of masses for which the scalar is stable in the 
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asymptotic region but unstable in the near- horizon region,^ For a given in this 
range, solutions exist in the region T//i < 7(m^), where 7(m^) is a dimensionless 
number depending on (which we evaluate in the next section). The value of 7(m^) 
is remarkably small for all in the allowed range. For example, with = —6 (not 
particularly close to the limiting value = —5), we have 7 ~ .00006333. It would be 
interesting to understand better how this small dimensionless number emerges since the 
setup has no small parameters. From the bulk point of view it is presumably related 
to the warping between IR and UV regions of the geometryp] From the boundary 
viewpoint, the low critical temperature may be explained by the BKL scaling pT f l ^ fT^] 
near a quantum critical point ^ 

For a given T and /i, we can use ^ to evaluate the action for the solution and 
compare this with the action for the soliton and/or Reissner-Nordstrom solution with 
the same T and fi. We find that the action for the new solutions is always less than 
the action for the Reissner-Nordstrom solutions, and is also less than the action for 
the soliton solutions for chemical potential in a region /i > fic{T). Thus, the solutions 
with scalar field represent the equilibrium phase in the region T/ fi < •y, fi > fic(T), as 
shown in figures [T]- [3] above. 

The transition between the deconfined and superconducting phases is second order, 
while the transition between confined and superconducting phases is first order. The 
place where these phase boundaries meet represents a triple point for the phase diagram 
where the three phases (confined, deconfined, superconducting) can coexist. 



4.2 Critical temperature 

For fixed m^, the value of ijj{0) in the solutions increases from zero at T//i = 7, 
diverging as T//i — )• 0. Since ij) is small everywhere near T/fi = 'y, the critical value of 
T/n will be the value where the ip equation, linearized around the Reissner-Nordstrom 
background, has a solution with the correct boundary conditions. Thus, we consider 
the equation 

r + (^ + ^)^' - "^v- = , (19) 

r g g 

where (setting r_|_ = 1) 

and find the value ^ = fic for which the equation admits a solution with boundary 
conditions ipi^) = 1 (we are free to choose this), = m? / g'{l) and the right falloff 

Solutions of this form were first found in lower dimensions in [14]. The zero-temperature limit of 
such solutions were considered in [25'. 

""^^By considering the alternate quantization mentioned in section 2 and fine-tuning the mass so that 
the dual operator has the smallest possible dimension consistent with unitarity in the dual field theory, 
we can obtain 7 as large as 0.0151, so even under the most favorable circumstances, the critical T//i 
is quite small. 

^^We thank D.T. Son for pointing this out to us. 
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Figure 5: Critical T//i vs of neutral scalar (filled circles). Mass is above BF 
bound asymptotically but below BF bound in near-horizon region of zero-temperature 
background solution in the range —6.25 < < —5. Unfilled circles represent critical 
values in the theory with alternate quantization of the scalar field, possible in the range 
-6.25 < < -5.25. 

{i/ji = 0) at infinity^ 

The choice r_|_ = 1 implies that T = (5 — 9/i^/8)/(47r), so we have 7 = (5 — 
9/i^/8)/(47r/ic). The results for 7(771^) are plotted in figure [s} For comparison, we also 
considered the theory defined with the alternate quantization {ip^ = 0) of the bulk 
scalar field (mentioned in section 2). As we see in figure [5| the critical temperatures 
are somewhat larger in this case, but still much smaller than 1 relative to /j,. 

4.3 Properties of the superconducting phase 

In the superconducting phase, it is interesting to ask how the charge density and free 
energy behave as a function of chemical potential. Since the solutions (as for the 
planar RN-black hole solutions) are trivially related to solutions where the w direction 
is non-compact, and since the underlying theory has a conformal symmetry, physical 
quantities in this phase (or in the RN phase) behave as fi^F{T/fi) for some non-trivial 

^■^To obtain a very accurate result, we first find a series solution ipiow near r — 1 with ■0(1) = 1 
(we are free to choose this) and = rn? /g'{l) and find a series solution iphigh for large r with the 

correct fall-off {tpi = 0) at infinity. Starting with ^low and '4''iow some r = n where the low r series 
solution is still very accurate, we then numerically integrate up to r = r2 where the large r series is 
very accurate and then find ^ for which -0^„„(r2)/-0„„„i(r2) = ^Lgh/'^high- 
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function F and a power At the critical value of T//i, we have a second order 
transition from the RN phase to the phase with scalar, so the free energy and its 
derivatives, and other physical quantities such as the density, are continuous across the 
transition. Thus, the relevant function F in these cases will be the same for the two 
phases across the transition. We find that the function F for either the charge density 
or the free energy changes very little between the very small value of T//i where the 
the transition occurs and the T — )■ limit. Thus, to a good approximation, we find 
that the density and free energy behave in the superconducting phase in the same way 
as for the zero temperature limit of the RN phase. For R=2/5, we have 

p ^ 0.320/ , 

while 

G ^ -.06V • 

In both cases, the behavior at large fi is governed by the underlying 4+1 dimensional 
conformal field theory. 



5 Charged scalar field: fiavor superconductivity 

In this section, we generalize our holographic model to the case where the scalar field 
is charged under the gauge field in the bulk. As we discussed in the introduction, this 
implies that the dual field theory includes some low-dimension gauge-invariant operator 
with charge, so the charge in this case is more naturally thought of as some isospin- 
type charge (since the smallest gauge-invariant operators carrying baryon charge have 
dimensions of order N). 

A significant qualitative difference in this case is that a scalar field condensate acts 
as a source for the electric field in the bulk, so it is possible to have solutions with no 
horizon carrying a finite charge density in the field theory. This gives the possibility 
of a fourth phase in which the scalar field condenses in the soliton background. 

To obtain the action for the charged scalar case, we begin with the action ^ 
and make the replacement d^ip — )■ d^ip ~ iQ^iiip- The results of the previous section 
correspond to g = 0. 



5.1 Low-temperature horizon free solutions with scalar 

Above some critical value of /i, there exist horizon-free geometries with a scalar field 
condensate. The solutions may be parameterized by the magnitude of the scalar at 
the IR tip of the geometry, and we will find a single solution for each such value. To 
determine these geometries, we need to take into account back-reaction on the metric. 
The most general solution with the desired properties can be described by the ansatz 

ds^ = r\e^^''^B{r)dw'' + dx^ + dy^ + dz^ - e^^^^rft^) ^ ^ 



^ If the solutions instead depended on the circle direction in a non-trivial way, we might have a 
general function of RT and 
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At = 0(r) , 



(20) 



where we demand A{oo) = C(oo) = and -B(oo) = 1. As for the sohton geometry, we 
expect that the w circle is contractible in the bulk so that B{rQ) = for some tq. For 
the geometry to be smooth at this point, the periodicity of the w direction must be 
chosen so that 

27iR = , , , . (21) 

Starting from the action (|2| with scalar derivatives replaced by covariant derivatives, 
the scalar and Maxwell equations are: 

„ /6 A' B' C'\ , 1 (e-^(q4>f ,\ , , 



W 2 B 2 r r^B ^ ' 



Following [19], we find that the Einstein equations give: 

_ 2r'^C" + r^C'^ + ArC + Ar^f^ - 2e-^(j)'^ 
~ r(8 + rC") ' 

5 --^ -o^^ +^7^ + ^ +^ J + ^ V' - 20 =0. 



r 2 I \ 2 2r*^ / \ 



These equations have two scaling symmetries, 



(26) 



?/'(r) = ^{ar) , 0(r) = -(j)(ar) , y4(r) = 74(ar) , _B(r) = B[ar) , C('") = C(ar) 

(X 



(27) 



and 



C' = C + A, = et0. (2^ 



Numerical evaluation of solutions 

To find solutions, we first use the scaling symmetries to fix rg = 1 and C(ro) = 0. For 
each value of ^/'(l), we use as a shooting parameter, choosing the value so that ip 
has the desired behavior for large r. From the solution obtained in this way, we can 



use (28) with A = — C(oo) to obtain the desired boundary condition C(oo) = in the 



rescaled solution. From (21), we see that the choice ro = 1 corresponds to a periodicity 



for the w direction equal to 

^-R=^T^- (29) 
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Figure 6: Action vs chemical potential for soliton with scalar solutions, taking = —6 
and q = 2. 



which will generally be different for solutions corresponding to different values of 
In order to obtain solutions corresponding to our chosen value R = 2/5 (such that the 
action for the soliton solution is -1) we use the scaling (27), taking a = B'(l)/5e~'^^°^''^'^. 
After all the scalings, we calculate the chemical potential and action (making use of 
(§) as 

/i = 0(cx)), S=[B]j^. 

The action is plotted against chemical potential for various values of q in figures |6j 
[7| and [s] taking the example of a mass just above the BF bound, = —6. 

We find that for large enough values of q, the chemical potential increases mono- 
tonically and the action decreases monotonically as we increase ipi^o)- This implies 
that we have a second order transition to the superconducting phase at a critical value, 
which can be determined by a linearized analysis (see appendix A) to be /i ^ 1.0125/g. 

Below q ^ 1.35, the chemical potential is no longer monotonic in ip{ro). We see 
that for q = 1.3, this results in a second order phase transition at /i ~ 1.558, followed 
by a first order phase transition at ~ 1.616 (taking R = 2/5). For smaller q (e.g. 
q = 1.2 in figure [s]), we simply have a first order transition to the superconducting 
phase at a value of chemical potential that is less than the value for the solution with 
infinitesimal scalar field. All of these results are completely analogous to the lower- 
dimensional results of 
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Figure 7: Action vs chemical potential for soliton with scalar solutions, taking = —6 
and q = 1.3. 
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Figure 8: Action vs chemical potential for soliton with scalar solutions, taking = —6 
and q = 1.2. 
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5.2 Hairy black hole solutions 

At high temperatures, the w circle is no longer contractible, and we assume that (as for 
the solutions without scalar field) the solution can be obtained by periodic identification 
of a solution with boundary i?^'^ instead of R^'^ x S^. Thus, we take the ansatz 

ds^ = -g(r)e-^('Ut^ + ^ + r\dw^ + dx^ + dy^ + dz^) , 

9{r) 

A = 0(r) , 
ip = ip{r) . 

The scalar and Maxwell's equations are 

,"+(i-| + ^),/ + i(^-™=).<^0. (30) 
V ^ g ) 9 \ 9 J 

^"^(^^)^'-^^-o. (31, 

while the Einstein equations are satisfied if 

j'+0-|)j+^ + ^-5r = O. (33) 

These have two symmetries: 

^(r) = ^(ar) , 0(r) = -0(ar) , xi^) = x{ar) , g{r) = ^g^ar) , (34) 

a 



A 



and 

X = X + A, = e-^0. (35) 

As we did for g = 0, we would like to find solutions with a horizon at some r = r+. 
The electric potential must also vanish at the horizon, and we are looking for solutions 
for which the leading falloff ipi in ^ vanishes for the scalar. Also, multiplying the first 
equation (30) by g and evaluating at r = r+ fixes ip'^r^) in terms of ipir^) and g'{r^). 
Altogether, our boundary conditions are 

g{r+) = 0, 0(r+) = O, X(oo)=0, = , 

and 



40r+ - 2m2r^V^(^+) - r+ex(^+)(0'(r+))2 ' 
The remaining freedom to choose r+ and 0'(r_|_) leads to a family of solutions with 
different T and /j,. Explicitly, we have 

= 0(oo) , T = ^g'ir+)e-^^^+y' . 



Solutions with the same T//x are simply related by the scaling symmetry (34). 
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Numerical evaluation of solutions 



To find solutions in practice, we can make use of the symmetries (34, 35) to initially 
set r_|_ = 1 and x(0) 



9(1) =0 



= and solve the equations with boundary conditions 
X(0) = 0, 0(1) = 0, <P'{l) = Eo, m=^o 



and 



F2 



0, and find one solution for each 



40 - 2m?^l 

We use Eq as a shooting parameter to enforce ipi 
ipQ. From these solutions, we apply the symmetry (35) with A = — x(oo) to restore 
x(cxd) = and finally use the symmetry (34) to scale to the desired temperature or 
chemical potential. 



5.3 Phase diagrams 

At a generic point in the phase diagram, we can have up to four solutions (AdS soliton, 
planar RN black hole, soliton with scalar, black hole with scalar), or more in cases 
where there is more than one solution of a given type. 

To map out the phase diagram, we evaluate the action for the various solutions 
using the methods of section 2. The equilibrium phase corresponds to the solution 
with lowest action. The phase diagrams for q = 1.3 and g = 2 (in the case = —6) 



are shown in figures 9 and 10/11) 



For large q, the condensation of the scalar field occurs in a region of the phase 
diagram where the back-reaction is negligible, so the phase diagram may be understood 
here (for /x ~ 1/g) by treating the gauge field and scalar on a fixed background (the 



Schwarzschild black hole). The resulting phase diagram is shown in figure (12). 



6 Discussion 

In this note, we have investigated the phase structure for a simple class of holographic 
systems which we have argued have the minimal set of ingredients to holographically 
describe the phenomenon of color superconductivity. Even in these simple models, we 
find a rich phase structure with features similar to the conjectured behavior of QCD 
at finite temperature and baryon chemical potential. It would be useful to verify the 
thermodynamic stability (and also the stability towards gravitational perturbations) 
of the phases that we have identified. This could indicate regions of the phase diagram 
where we have not yet identified the true equilibrium phase for the model, for example 
since our ansatz might be too symmetric. 

We have calculated some of the basic thermodynamic observables, but it would be 
interesting to investigate more fully the physical properties of the various phases and 
establish more definitively a connection between the phase we find at large /i and small 
temperature and the physics of color superconductivity. 
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Figure 9: Phase diagram for = —6 and q = 2. Clockwise from the origin, the 
phases correspond to the AdS sohton (confined), RN black hole, black hole with scalar, 
and soliton with scalar. 




Figure 10: Phase diagram for = —6 and q = 1.3. Clockwise from the origin, the 
phases correspond to the AdS soliton (confined), RN black hole, black hole with scalar, 
and soliton with scalar. 
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Figure 11: Small temperature region of phase diagram for = —6 and q = 1.3. 
Dashed line represents a first order transition within the soliton with scalar phase. 




Figure 12: Phase diagram for large q, m? = —6. 
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Apart from the ipipip'^ip'^ condensate that we can see directly using the ingredients of 
our model, there are various other features that characterize a color superconductivity 
phase [I] . Typically, the breaking of gauge symmetry is accompanied by some breaking 
of exact or approximate flavor symmetries. Thus, the superconducting phase has a 
low-energy spectrum characterized by Goldstone bosons or pseudo-Goldstone bosons 
associated with the broken flavor symmetries, together with massive vector bosons 
associated with the spontaneously broken gauge symmetry. It would therefore be 
interesting to analyze the spectrum of fluctuations in our model to compare with these 
expectations. 

A caveat related to looking for features associated with the global flavor symmetries 
(and their breaking) in our model is that we may not have included enough ingredients 
in our bottom-up approach for all these features to be present. In simple models where 
the flavor degrees of freedom are associated with probe branes, there are explicit gauge 
fields in the bulk dual to the global symmetry current operators. However, in fully back- 
reacted solutions (appropriate for studying Nj ~ Nc), these branes are replaced by a 
modified geometry with additional fluxes (for an explicit example of such solutions, see 
|29j). In these solutions (which we are trying to model in our approach), it is less clear 
how to identify the global symmetry group from the gravity solution, but presumably 
it has to do with some detailed properties of the geometry. Thus, it is possible that the 
Goldstone modes associated with broken flavor symmetries correspond to fluctuations 
in some fields (e.g. form-fields) that we have not included. 

The color superconducting condensate also breaks the global baryon number sym- 
metry, so there should be an associated Goldstone boson related to the phase of the 
condensate, and associated superfluidity phenomena. In other holographic models 
with superfluidity, the condensate is dual to a charged scalar field in the bulk and the 
Goldstone mode is related to fluctuations in the phase of this field. However, as we 
mentioned in the introduction, the baryon operator has dimension of order A^, so we do 
not expect a light charged scalar field in the bulk. In a more complete top-down model, 
the baryon operator may be related to some non-perturbative degrees of freedom (such 
as D-branes) in the bulk, and it may be necessary to have a model with these degrees 
of freedom included in order to directly see the Goldstone mode from the bulk physics. 
Related to these observations, it may be interesting to probe our model with D-branes 
(put in by hand), in order to make the relation to microscopic physics more manifest, 
and to help gain a better understanding of the phenomenological parameters of our 
model. 

There are a number of variants on the model that would be interesting to study. 
First, the breaking of scale-invariance, implemented in our model by the varying cir- 
cle direction in the bulk, could be achieved in other ways, replacing (7^^ with a more 
general scalar field, as in the model of |Tl]. In the setup of that paper, the transi- 
tion between confined and deconfined phases was found to exhibit crossover behavior 
at small chemical potential, a feature expected in the real QCD phase diagram and 
expected generally for massive quarks with sufficiently large Nf/Nc- It would be inter- 
esting to look for an even more realistic holographic model by incorporating features 
of the model we have studied here and the model of 
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It would also be interesting to look at the effects of a Chern- Simons term for the 
bulk gauge field. In [21] and [22], it was shown that such a term (with sufficiently 
large coefficient) gives rise to an instability toward inhomogeneous phases, perhaps 
associated with the chiral density wave phase believed to exist at large density in 
QCD with iVc <^ A^/ [El [7|. It is interesting to investigate the interplay between 
these inhomogeneous instabilities and the superconducting instabilities discussed in 
the present paper. It would also be interesting to consider more general actions (such 
as Born-Infeld) for the gauge field, interaction terms for the scalar field in the bulk, or 
other couplings between the scalar field and gauge field. 

Finally, once the technical challenges of writing down fully back-reacted solutions 
for top-down models of holographic QCD with Nf Nc have been overcome, it will be 
interesting to see whether the basic features we find here are manifested in the more 
complete string-theoretic models. If certain features are found to be universal, these 
might taken as qualitative predictions for the QCD phase diagram, or at least motivate 
an effort to understand whether these features are also present in the phase diagram 
of real- world QCD. 
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A Large charge limit 

In this appendix, we analyze the case of large q. This is particularly simple, since in 
this limit, the back-reaction of the scalar and the gauge field on the metric go to zero in 
the region of the phase diagram where transitions to the superconducting phases occur. 
Explicitly, we can show that in the limit g — t- oo with qfi fixed, the gauge field and 
scalar field decouple from the equations for the metric, but still give rise to a nontrivial 
phase structure. To investigate this, we need only consider the scalar field and gauge 
field equations on the fixed background spacetimes corresponding to low temperatures 
(the soliton geometry) and high temperatures (the Schwarzschild black hole). 

Low Temperature 

Starting from the action ^ for the scalar field and gauge field on the soliton background 
([7|, we find that the equations of motion are (setting L = 1) 
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where / is defined in ([8j). 

These equations have two scahng symmetries related to the conformal symmetry 
of the boundary field theory and to the absence of back-reaction in our large charge 
limit. Given a solution (0(r), '?/'(r), tq, g, m), we can check that the scaling 

(0(r), V^(r), ro, g, m) (/30(ar), (^aipiar), — , m) 

a ap 

sends solutions to solutions. For our calculations, we will use this to set rg = g = 1. 

Multiplying these equations by / and taking the limit r — )■ ro = 1, we find that 
regular solutions must obey 



0'(1) 



2V^2(i)0(i) 



5 

We have two remaining parameters, V^(0) and 0(0). One of these can be fixed by 
demanding that the "non-normalizible" mode of vanishes at infinity, while different 
values of the remaining parameter correspond to different values of /i. 

Employing numerics, we find that for a fixed value of m^, there is some critical 
value of fi above which solutions with a condensed scalar field exist. 

In order to determine the critical value /ic('^^), we use the fact that the field values 
go to zero as we approach the critical /i from above. Thus, at the critical /i, the 
equations above linearized around the background solution = /i should admit a 
solution with the correct boundary conditions. The linearized equations decouple from 
each other, so we need only study the ip equation. This becomes 

-ijj" ( ^^^'^ ~ "^^'"^)^ g 

\r{r^ — I) J — 1 

We can take '0(1) = 1 without loss of generality, so the boundary condition for ip' 
becomes 

^'(1) = J(m2-/i2). 

Given m^, we now find /x^ by demanding that the leading asymptotic mode (0^1) of 



vanishes. Our results for the critical as a function of m are shown in figure 13 



High temperature 

The high temperature geometry relevant to the limit of large q with fiq fixed is the 
/i — )■ limit of the Reissner- Nordstrom geometry (10), which gives the planar AdS- 
Schwarzschild black hole (with one of the spatial directions compactified) . This is the 
relevant background for T > l/(27ri?). 
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Explicitly, we have 

ds'^ = (-dt^fir) + dx^ + dy^ + dz^ + dw 



dr 



2 



where 

Here, r+ is related to the temperature by 

AttT 
'^ = ^- 

The equations of motion in this background are 

The equations have the same scaling symmetry as before, so we can set r+ = g = 1 for 
numerics. Here, the choice r+ = 1 corresponds to T = l/{27iR), where R is the radius 
chosen in the previous section by setting tq = 1. In this case, the boundary conditions 
are 

0(1) = , ^'{1) = ^ . 

To determine the physics at other temperatures, we can fix q and R and use the scaling 
to adjust the temperature. 

For any values of parameters, we have a solution 



^ = , 0(r) = 



r 



3' 



corresponding to the pure Reissner-Nordstrom background in the probe limit. 

As in the low temperature phase, we find a critical value /ic = F{m?) (or, restoring 
temperature dependence, Hc = ^F{m'^)) for each choice of m^, above which there 
is another solution with nonzero ip. This critical fi may again be determined by a 
linearized analysis, from which we obtain the equation 

We can set il^iX) = 1 without loss of generality, and this requires 

^.'(1) 3 . 

These can be solved numerically to find F{m?)^ and our results (with the low temper- 



2 

m 



ature results) are plotted in figure 13 



A sample phase diagram, for the case m = —6 is shown in figure 12 
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Figure 13: Critical values of jj^q vs for scalar condensation in large q limit. The 
top curve is the critical value for /i in black hole phase (just above the transition 
temperature), while the bottom curve is the critical in low temperature phase. 

A.l Order of phase transitions in the probe Umit 

To complete this section, we verify analytically that the action for solutions with scalar 
field in the probe limit is always less than the corresponding unperturbed solution. In 
this limit we neglect the gravity back reaction of the gauge fields and scalar. The 
on-shell action in this approximation is given by 

^ = J d'^'x V^gg'V^^ . (36) 

We have used the fact that the scalar action is quadratic and vanishes on-shell once 
the boundary value of scalar is kept to zero |24j . Writing the action in this simple form 
gives us information about the relative free energy of the different phases. 
The solution for At in the superconducting phase may be written as 

Af = A° + 6 At , (37) 

where 6 At — )■ in the IR region of the bulk and near the boundary. A^ is the value of 



At in the normal phase. Then, from eq. (36) we get 

= ^ + ^ j drV^99'^9'%A%{5At) + j V^gg-g^'^^dr . (38) 

The cross term between A^ and 5 At vanishes after integrating by parts and then using 
the eom of A?. Hence 



5S = - Sold = (TV) j ^gg"g''^^^dr < Oas^?" < 0. (39) 

Therefore if a phase with non-trivial scalar condensate exists it will always have a lower 
free energy than the normal phase and the associated transition will be of second order. 

The introduction of gravity may give rise to a positive term in the on-shell action 
and the nature of phase transition may change. 
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Figure 14: Critical T//x vs charge q for condensation of m? 
Reissner-Nordstrom background. 



-6 scalar field in 



B Critical /i for solutions with infinitesimal charged 
scalar 

To find the critical at which solutions with infinitesimal scalar field exist, we find the 
value of /i for which the linearized scalar equation about the appropriate background 
admits a solution with the right boundary conditions at infinity. 
At low temperatures, this gives (setting tq = 1) 

g(r) = - 4j , = /i , 



while for the RN black hole background (setting r 

g{r) = - (^1 + 5 
= ^ (l - 



1) we have 



8r6 ' 



More general values of Tq or r+ can be restored by the scaling symmetry. 

For = —6, we find a critical value of in the low-temperature case given by 
l^iowQ. = 5.089/(27ri?). At high temperatures, the critical solutions exist T/fi when has 



a critical value as plotted in figure 14 
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